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ABSTRACT 

Recent data on the high-redshift abundance of damped Lya systems are compared with 
theoretical predictions for 'blue' (i.e. n > 1) Mixed Dark Matter models. The results 
show that decreasing the hot component fraction n v and/or increasing the primordial 
spectral index n implies an earlier epoch of cosmic structure formation. However, we 
also show that varying Q v and n in these directions makes the models barely consistent 
with the observed abundance of galaxy clusters. Therefore, requiring at the same time 
observational constraints on damped Lya systems and cluster abundance to be satisfied 
represents a challenge for the Mixed Dark Matter class of models. 



1 INTRODUCTION 

Since a long time observations of high-redshift objects have 
become a potentially powerful constraint for models of cos- 
mic structure formation. The availability of statistically re- 
liable samples of quasars allowed to address this problem in 
a quantitative way in the framework of the Cold Dark Mat- 
ter cosmogony (Efstathiou & Rees 1988; Haehnelt 1993). 
Moreover, the comparison of predictions and observations of 
quasar abundance at different redshifts has been used as a 
test for model reliability (e.g. Nusser & Silk 1993; Pogosyan 
& Starobinsky 1993). 

Recently, damped Lya systems (DLAS) have been rec- 
ognized as a promising way to trace the presence of high 
redshift collapsed structures, thanks to the possibility of 
identifying them as protogalaxies and to their detectabil- 
ity at high z (see Wolfe 1993 for a comprehensive review). 
DLAS are seen as wide absorbtion features in quasar spec- 
tra. The associated absorbing systems have a neutral hydro- 
gen column density > 10 20 cm -2 . The rather large abun- 
dance of DLAS makes it possible to compile reliable statis- 
tical samples (Lanzetta 1993; Lanzetta, Wolfe & Turnshek 
1995; Storrie-Lombardi et al. 1995; Wolfe et al. 1995). Once 
the parameters of the Friedmann background are specified, 
observations of DLAS can be turned into the value of the 
cosmological density parameter Q g contributed by the neu- 
tral gas, which is associated with DLAS. It turns out that at 
z ~ 3 this quantity is comparable to the mass density of vis- 
ible matter in nearby galaxies, thus suggesting that DLAS 
trace a population of galaxy progenitors. 

Based on the APM QSO catalogue, Storrie-Lombardi 
et al. (1995) recently presented the most extended DLAS 
sample up-to-date, covering the range 2.8 < z < 4.4. In the 
following of this paper we will consider their highest redshift 



data as the most constraining ones and we will compare 
them with model predictions. 

Several authors (Subramanian & Padmanabhan 1994; 
Mo & Miralda-Escude 1994; Kauffmann & Chariot 1994; 
Ma & Bertschinger 1994) have recently claimed that the 
large value of Q g observed at z ^ 3 is incompatible with pre- 
dictions of the Mixed (i.e. cold+hot) Dark Matter (MDM) 
model with spectral index n = 1 and Q. v ~ 0.3 contributed 
by one species of massive neutrinos and fij, =0.1 for the 
baryon fractional density (Klypin et al. 1993; Nolthenius, 
Klypin & Primack 1995). Klypin et al. (1995) reached sub- 
stantially the same conclusions about this model, but em- 
phasized two relevant points: (?) any theoretical prediction 
is very sensitive to the choice of the parameters of the 
model needed to obtain il g from a given power-spectrum; 
(ii) slightly lowering fi„ to 20-25% keeps MDM into bet- 
ter agreement with DLAS data, independently of whether 
the hot component is given by one or two massive neutrino 
species (see also Primack et al. 1995). 

A possible alternative, still in the framework of MDM, 
is to consider 'blue' (n > 1) primordial spectra of density 
fluctuations. The advantage of these models is an anticipa- 
tion of the epoch of structure formation due to the higher 
small-scale power. The choice of blue spectra was originally 
suggested by the analysis of Cosmic Microwave Background 
anisotropies on scales larger than 1° (e.g. Devlin et al. 1994; 
Hancock et al. 1994; Bennett et al. 1994). Possible indica- 
tions for blue spectra comes from large bulk flows (Lauer & 
Postman 1994; see however Riess, Press & Kirshner 1995; 
Branchini & Plionis 1995) and large voids in the galaxy dis- 
tribution (Piran et al. 1993). In recent years many authors 
have pointed out that the inflationary dynamics can easily 
account for the origin of blue perturbation spectra (Liddle 
& Lyth 1993; Linde 1994; Mollerach, Matarrese & Lucchin 
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Figure 1. The r.m.s. fluctuation amplitude within a top-hat sphere of Sh^ 1 Mpc, <rg, for COBE normalization (see text). Left panel: 
dependence on the hot component fraction Q v ; different lines refer to n = 0.9,1,1.1,1.2,1.3,1.4 from bottom to top. Right panel: 
dependence on the spectral index n; different lines refer to H„ = 0, 0.1, 0.2, 0.3, 0.4, 0.5 from top to bottom. 



1994; Copeland et al. 1994), in particular in the framework of 
the so-called hybrid models. Recently Lucchin et al. (1995), 
using linear theory and N-body simulations, performed an 
extended analysis of the large-scale structure arising from 
blue MDM (BMDM) models: the most interesting advantage 
of these models is the increase of the galaxy formation red- 
shift: for instance, taking Q. v = 0.3 one has for the redshift 
of non-linearity on galactic scale (M = 10 12 M©) z n i ~ 1.9 
if n = 1.2 and z ni m 0.6 if n — 1. The same class of mod- 
els has been tested against observational data, using linear 
theory predictions, by Dvali, Shan & Schaefer (1994) and 
Pogosyan & Starobinsky (1995a). 

In this work we will compare BMDM model predictions, 
for different values of 0„ and n (for the sake of comparison 
we also consider the 0.9 < n < 1 tilted models), with the ob- 
served DLAS abundance. Furthermore, we will also discuss 
the implications of the observed abundance of galaxy clus- 
ters for BMDM models. In fact, increasing n for a fixed Q. v 
and fixed normalization to COBE rises up the value of as, 
the r.m.s. fluctuation within a top-hat sphere of 8 h^ 1 Mpc 
radius, which is constrained by the observed abundance of 
galaxy clusters to be as ~ 0.6 (White, Efstathiou & Frenk 
1993). 



2 THE METHOD 

In order to connect our model predictions to DLAS observ- 
ables, let us define O co ;f(z) as the fractional matter density 
within collapsed structures at redshift z. Therefore 

GcouM = ^ , (l) 

where fl b is the fractional baryon density (since h = 0.5 is 
assumed throughout the paper, we take fij, = 0.05 accord- 
ing to standard primordial nucleosynthesis; see, e.g., Reeves 
1994) and f g is the fraction of the HI gas, which is involved 
in DLAS. Although the observed decrease of Q g with red- 
shift for 2 < 3.5 is usually considered as an indication of 



gas consumption into stars (e.g. Lanzetta et al. 1995; Wolfe 
et al. 1995), the actual value of f g at the high redshift we 
are interested in is not clear. In any case, since Q g at such 
high redshifts is quite similar to the fractional density con- 
tributed by visible matter in present-day normal galaxies, 
we expect f g not to be a particularly small number. In the 
following we will show results based both on f g = 0.5 and 
1. 

Taking h — 0.5 and an Einstein-de Sitter universe, the 
data at z = 4.25 from Storrie-Lombardi et al. (1995) turn 
into n co u = (8.8 ± 2.0) x 10~ 2 and (4.4 ± 1.0) x 10~ 2 for 
f g = 0.5 and 1, respectively. 

From the theoretical side, the Press & Schechter (1974) 
approach gives a recipe to compute the contribution to the 
cosmic density due to the matter within collapsed structures 
of mass M at redshift z: 

n co u(M,z) = erfcf Sc — V (2) 

The above expression assumes Gaussian fluctuations and 8 C 
is the linearly extrapolated density contrast for the collapse 
of a perturbation; om is the r.m.s. fluctuation at the mass- 
scale M, where 

M = (2tyR 2 ) 3/2 P (3) 

for the Gaussian window that we will assume in the follow- 
ing. Here, p is the average matter density, which is taken to 
have the critical value. 

As for the mass of the structures hosting DLAS, it has 
been argued that, since the high column density of the ab- 
sorber is typical of large disks of luminous galaxies, DLAS 
should be located within massive structures of ~ 10 12 M . 
However, it is not clear at all whether the properties of 
present-day galaxies can be extrapolated to their high- 
redshift progenitors. Therefore, we prefer to leave open the 
possibility that DLAS are hosted within smaller structures. 
It is clear that, when a model is in trouble in accounting 
for the DLAS abundance if the hosting structure is a dwarf 
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Figure 2. The fractional matter density within collapsed structures £l co ii at redshift z = 4.25, when 8 C = 1.5 and M = 10 11 Mq are 
assumed. Left panel: dependence on the hot component fraction £7„; different lines refer to n = 0.9, 1, 1.1, 1.2, 1.3, 1.4 from bottom to 
top. Right panel: dependence on the spectral index n; different lines refer to Q v = 0, 0.1, 0.2, 0.3, 0.4, 0.5 from top to bottom. Error bars 
show the effect of taking M = 10 10 Mq and 10 12 Mq. The horizontal lines refer to the observational data by Storrie-Lombardi et al. 
(1995) with the corresponding uncertainties. 



galaxy (M ~ 10 10 Mq), it would certainly be ruled out if 
more massive protogalaxies are required. 

Linear theory for the top-hat spherical collapse pre- 
dicts S c = 1.69. However, effects of non-linearity as well 
as asphericity of the collapse could cause significant devia- 
tions from this value. Klypin et al. (1995) estimated the halo 
abundance at different redshifts from high mass resolution 
N-body simulations. By using the Gaussian window, they 
found a good agreement with the Press-Schechter expres- 
sion for values as low as S c = 1.3-1.4 (see also Efstathiou & 
Rees 1988). On the other hand, Ma & Bertschinger (1994) 
found that for Q. v = 0.3 and a top-hat window S c ~ 1.8 
is always required, which corresponds to S c ~ 1.7 for the 
Gaussian window. In the following we prefer to show results 
based on S c = 1.5 but, due to the previous uncertainties, we 
will discuss also the effect of different choices for 5 C . 

For the MDM transfer function we take the fit obtained 
by Pogosyan & Starobinsky (1995a), which provides a con- 
tinuous dependence on the fractional density Q. v contributed 
by one massive neutrino. As for the cold part of the trans- 
fer function, we use the Cold Dark Matter expression by 
Efstathiou, Bond & White (1992), with the shape param- 
eter r = Q h exp(— 2f2;,), according to the prescription of 
Peacock & Dodds (1994), to account for the baryonic com- 
ponent. We varied fl„ in the interval < Q v < 0.5. We as- 
sume for the primordial (post inflationary) power-spectrum 
P(k) oc k n with 0.9 < n < 1.4. Each model is normalized to 
the 9-th multipole component of the COBE DMR two-year 
data, 09 = 8.2, which has been shown to be independent of 
the n value to a good accuracy (Gorski et al. 1994). 

In Figure [l] we plot the resulting erg value. As for the f2„ 
dependence, results are plotted in the left panel for n = 0.9— 
1.4 going from lower to upper curves with steps of 0.1. In a 
similar fashion, in the right panel we plot the n dependence. 
Going from higher to lower curves, we plot results for Q. v = 



— 0.5 with steps of 0.1. As expected, as is an increasing 
function of n, while it decreases with Q, v . 



3 DISCUSSION 

The results of our analysis on DLAS are summarized in 
Figure ^, where we plot tl co u, estimated at z = 4.25, as 
a function of Q v (left panel) and of n (right panel), after 
assuming 5 C = 1.5 and M — 10 11 Mq. In each panel, dif- 
ferent curves are for the same choice of parameters as in 
Figure [1} Upper and lower error bars show the effect of tak- 
ing M = 10 10 M and 10 12 Mq, respectively. The horizontal 
solid line is the observational result with the corresponding 
uncertainties (dotted lines) , which is obtained by converting 
the Q g value, as reported by Storrie-Lombardi et al. (1995) 
at z — 4.25, to Q, co ii according to eq.([l|) with f g = 1. 

Figure ^ shows in the n u - n plane the models which 
reproduce the observed Cl C oll> taking f g = l (left panel) and 
f g = 0.5 (right panel). The heavy solid curve corresponds 
to 5 C = 1.5 and M = 10 11 Mq, with the lighter curves 
delimiting the observational uncertainties. Upper and lower 
dashed lines show the effect of varying S c to 1.3 and 1.7, 
respectively. Upper and lower dotted curves refer to M = 
10 10 and 10 12 Mq, respectively. The overall result that we 
get is that decreasing the hot component fraction Q. v and/or 
increasing the primordial spectral index n implies an earlier 
formation of cosmic structures. 

It should be noted that realistic observational uncer- 
tainties should be larger than the error bars reported by 
Storrie-Lombardi et al. (1995), since they do not include 
any systematic observational bias. Recently, Bartelman & 
Loeb (1995) emphasized the role of the amplification bias, 
due to DLAS gravitational lensing of QSOs, in the DLAS 
detection. They pointed out that (a) lensing effects bias 
upwards the 0. g {z) value by an amount depending on the 
parameters of the Friedmann background, as well as on the 
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Figure 3. The models in the £7„ — n plane reproducing the observed £l co u, taking f g = 1 (left panel) and f g = 0.5 (right panel). The 
heavy solid curve is for S c = 1.5 and M = 10 Mq; lighter solid curves refer to the observational uncertainties. Dashed and dotted lines 
show the effect of varying 8 C and M, respectively (see the text for the assumed values). 



redshift; (b) the observed absorber sample may be biased 
toward larger values of their internal line-of-sight velocity 
dispersions leading to an overestimate of the total absorber 
mass. Therefore, both effects go in the direction of allevi- 
ating the galaxy formation redshift problem. On the other 
hand, Fall & Pei (1995) detailed the consequences of dust ab- 
sorbtion in DLAS. They argued that dust obscuration causes 
incompleteness in the optically selected quasar samples, and, 
therefore, in the DLAS samples as well. In this case, the re- 
sulting Q g (z) is biased downwards by an amount depending 
on the model for the DLAS chemical evolution. 

It is however clear that, even taking the observational 
results at face value with their small error bars, the rather 
poor knowledge of the parameters entering in the Press- 
Schechter prediction for Q g (i.e. S c , M and f g ) makes it 
difficult to put stringent constraints on fi„ and n. 

For instance, if one takes 1.3 ;$ 5 C ^ 1.5, as suggested 
by several N-body simulations (e.g. Efstathiou & Rees 1988; 
Klypin et al. 1995) and analytical considerations on the 
Press-Schechter approach (e.g. Jain & Bertschinger 1994), 



0.2 and n = 1 would be allowed for M 



in 1 



10 Mq, unless f g is sensibly below unity. On the other 
hand, blueing the spectrum to n = 1.2 increases the allowed 
hot fraction to Q. v > 0.4, unless S c ~ 1.7 or M ~ 10 12 M 
are taken. In order to more tightly constrain the models, a 
better understanding of galaxy formation through hydrody- 
namical simulations would be needed to clarify what DLAS 
actually are. This would provide more reliable values for S c , 
M and f g . 

It is however clear that those models which fit the data 
at z ~ 4 need also to be tested against present-day observ- 
ables. One of these tests is represented by the abundance of 
galaxy clusters, which has been shown to represent a pow- 
erful constraint for dark matter models (e.g., White et al. 
1993). In the Press & Schechter (1974) approach, the num- 
ber of density of clusters with mass above M is given by 



where 
n(M) dM 



exp 



dR 



2a 2 {R) I R? 



(5) 



JV(> M) = I n(M')dM' , 

' M 



(4) 



is the average cluster number density with mass in the range 
[M, M + dM] . In the above expression, the quantities 

a 2 {R) = ^_ J k 2 P{k)W 2 {kR)dk , 

n(R) - 1 [ ^P(k) dW2{kR) dk (6) 

V{H) ~ 2>kW{R) J k F{k) d(kR) kR [b > 

convey the information about the power-spectrum. As be- 
fore, we use a Gaussian window for W(kR), so that a — 
(27r) 3 / 2 in eq.(^) and the mass M is related to the scale 
R according to eq.(^|). Klypin & Rhee (1994) found that 
8 C ~ 1.5 for their MDM cluster N-body simulations, with 
Q v — 0.3 and Gaussian filter (see Borgani et al. 1995, for the 
dependence of the cluster mass function on S c for different 
dark matter models). 

As for observational data, White et al. (1993) esti- 
mated a cluster abundance of about 5 x 10~ 7 Mpc -3 for 
masses exceeding M = 8.4 x 10 14 Mq using X-ray data. 
Biviano et al. (1993) based their analysis on observed clus- 
ter velocity dispersion and obtained an abundance of about 
7.5 x 10 -7 Mpc -3 for clusters exceeding the above mass 
limit. 

In Figure ^ we plot the fi„ - n relations for different val- 
ues of N(> M) (see caption) and taking M = 8.4 x 10 14 M . 
Quite remarkably, for S c = 1.5 (left panel) no values of tt v 
and n in the considered ranges give a cluster abundance 
as low as the observational ones. For instance, assuming 
(£2„,n) = (0.2, 1) it is N(> M) ~ 3.4 x 10 -6 Mpc" 3 , while 
N(> M) ~ 7.1xl0~ 6 Mpc" 3 for (fi„,n) = (0.3, 1.2). In gen- 
eral, lowering Q. v and/or increasing n, as suggested by the 
DLAS analysis, makes the disagreement even worse. This re- 
sult agrees with the expectation that as — 0.6 (White et al. 
1993) to reproduce the observed cluster abundance, while in 
general larger normalizations are required by our models (cf. 
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Figure 4. The models in the f2„ - n plane corresponding to different values for the cluster abundance, N(> M), taking M = 8.4 X 
10 14 Mq. Left and right panels correspond to assuming S c = 1.5 and S c = 1.7, respectively. Solid, dotted, short— dashed, long-dashed and 
dot-dashed curves (partly absent in the left panel) are for log7V(> M) = —5.25, —5.5, —5.75, —6., —6.25, respectively. For comparison, 
the observational results by White et al. (1993) and by Biviano et al. (1993) give logiV(> M) ~ —6.3 and —6.1, respectively. 



Figure |l|). Increasing S c to 1.7 alleviates to some extent the 
disagreement. However, even in this case, in order to repro- 
duce the observational N(> M) one should have 0.2 
and 0.95, with the limiting values (f2„,n) = (0.2,0.9) 
being only marginally consistent with the DLAS constraints 
corresponding to most optimistic choice of the parameters 
(cf. the upper dotted curve in the left panel of Figure 2). 
This general problem of the models in reproducing the clus- 
ter abundance was also recognized by Pogosyan & Starobin- 
sky (1995a). The point is also discussed by Lucchin et al. 
(1995). 

Although systematic observational uncertainties could 
well affect the determination of cluster masses from both 
X-ray and velocity dispersion data, it is not clear whether 
they can justify the order of magnitude (or even more) dis- 
crepancy between the data and those models which would 
have been preferred on the ground of DLAS constraints. 

A first possibility to alleviate this problem would be to 
increase the baryon fraction Qi,. On one hand, this has the 
effect of lowering the small-scale fluctuation amplitude and, 
therefore, erg. On the other hand, this fluctuation suppres- 
sion should not damage too seriously DLAS predictions on 
Q g , since this effect is partly compensated by a larger de- 
nominator in eq.(jl|). For instance, taking 5 C — 1.5 for the 
(p,„,n) = (0.2,0.9) model the cluster abundance changes 
from n(> M) ~ 1.4 x 10" 6 Mpc -3 to n(> M) ~ 8.2 x 10" 7 
Mpc -3 when passing from 5% to 10% of baryonic fraction. 
However, such an effect turns out not to be effective in rec- 
onciling with observational data those models which largely 
overproduce clusters. Indeed, even for S c — 1.7, n(> M) for 
(fi„,n) = (0.3,1.2) drops only from 7.1 x 10~ 6 Mpc -3 to 
2.5 x 10 -6 Mpc -3 when passing from 5% to 20% of baryonic 
fraction, the second value already being largely inconsistent 
with the primordial nucleosynthesis predictions. 

A further possibility is sharing the hot component be- 
tween more than one massive neutrino species (Primack et 
al. 1995; Pogosyan & Starobinsky 1995b; Babu, Schaefer & 



Shaft 1995). The subsequent variation of the neutrino free- 
streaming has been shown to decrease as to an adequate 
level, without significantly affecting results at the galactic 
scale, which is relevant for DLAS. 

As a general conclusion we would stress the effective- 
ness of putting together different kinds of observational con- 
straints to restrict the range of allowed models. As we have 
shown, the effect of blueing the primordial perturbation 
spectrum goes in the direction of increasing the redshift of 
structure formation. However, this also increases the r.m.s. 
fluctuation on the cluster mass scale to a dangerous level. 
Deciding whether the narrowing of the allowed region of the 
fl u — n plane points towards the selection of the best model 
or towards ruling out the entire class of models requires a 
clarification of both the observational situation and of how 
models have to be compared to data. 
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